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The virial theorem for non-relativistic complex fields in D spatial dimensions and with arbitrary 
many-body potential is derived, using path-integral methods and scaling arguments recently de¬ 
veloped to analyze quantum anomalies in low-dimensional systems. The potential appearance of a 
Jacobian J due to a change of variables in the path-integral expression for the partition function 
of the system is pointed out, although in order to make contact with the literature most of the 
analysis deals with the J = 1 case. The virial theorem is recast into a form that displays the effect 
of microscopic scales on the thermodynamics of the system. From the point of view of this paper 
the case usually considered, J = 1, is not natural, and the generalization to the case J 7 ^ 1 is briefly 
presented. 
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I. INTRODUCTION 

The virial theorem has been proven using a variety of methods. Recently, a path-integral derivation of the virial 
theorem has been developed in the context of quantum anomalies in non-relativistic 2D systems, or more generally, 
systems with SO(2, 1) classical symmetry [I|. The path integral is most useful in isolating the anomaly contribution 
to the equation of state so obtained. This method is in fact quite general, and applicable for non-relativistic systems 
with an arbitrary 2-body potential V{xi,X 2 ) in D spatial dimensions, even when there are no quantum anomalies 
present. We present such derivation in this note, extending the original derivation using also diagrammatic analysis, 
and recasting the virial theorem into a general equation that relates macroscopic thermodynamics variables to the 
microscopic physics. As it will be shown, there is generically a Jacobian term J that may contribute to the virial 
theorem, regardless of the existence of a classical scaling symmetry. We will mainly concern ourselves here with the 
case J = 1 (which we term “non-anomalous”). Comments and conclusions end the note. 


II. VIRIAL THEOREM 

The work in [l| was based partly on the work by Toyoda et al. They postulated that spatial scalings^ 

X' = Xx , 

ifi'{t,x') = X~^^^'ip{t,x), 


leave the particle number density invariant: 


d^xip^t,x)ip{t,x) = d^x' d:' x'). 


( 2 ) 


*Electronic address: cllin@uh.edu Electronic address: cordonez@central.uh.edu 

^ Toyoda et al. introduced an auxiliary external potential that has the effect of confining the system to a volume V, and then, through a 
series of infinitesimal scalings and algebraic arguments derived what amounts to the equation of state, which they referred to as virial 
theorem. Unlike them, we’re not using an external potential but simply consider a system with a large volume U(so all the typical 
large-volume thermodynamical considerations apply), but like them, we’re also calling virial theorem the equation of state that will be 
derived in this paper. 
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Let us consider a non-relativistic system whose microscopic physics is represented by a generic 2-body interaction^ 

£ = V'* V'- ^ J d^yi’*{t,x)'ip{t,x)V{x-y)i}*{t,y)'tp{t,y). (3) 

Giving our system a macroscopic volume V, temperature /3“^, and chemical potential y, and going into imaginary 
time gives for the partition function: 


Z[V P fj]= J [dtlj*][dll}]e~ dT/v [b*(ax-^-M)b+5/d^yb*('r,S)b('r.^)V(£-y)j/'*('r,y)i/'('r,i/)] 

Now consider a new system with the same temperature and chemical potential, but at volume V = \^V'. 


Z[\^V 13 y] = J [d‘lp'*][dtp’]e~ 

Substituting Eq. (1) into Eq. (5) gives: 


Z[X^V,i3,y]= [d^*][d^]Je 


- So d-r Jv / d’^ y tjj* {T,x)tp{T,x)V{\(x -y) )tp* {T,y)ili(T,y)^ 


( 6 ) 


where J is the Jacobian for the transformation (■;/;'*,"0') —> (0*,0). As mentioned above, our emphasis will be 
in the non-anomalous case, and henceforth we assume J = 1 (see however comments and conclusions). Then 
Z[X^V, /3, y] = Z^[V, /3, y], where the superscript A represents a microscopic system whose kinetic energy has a factor 
^ and whose potential is V{X{x —y))- Note that Z^^^[V, P, y] = Z[V,P,y\. 


The pressures corresponding to Z[X^V, /3, y] and Z[V, /3, y] are equal, since the intensive variables y and are the 
same, and they correspond to the same microscopic system. The argument we just made for the pressures being the 
same is valid in the thermodynamic limit, based on the principle that two intensive variables determine the third via 
an equation of state e.g., P = pT for an ideal gas. However, in the next section we will also provide a diagrammatical 
proof that the two pressures are the same. 


For now assume the pressures are equal. Then using Z = , we get: 


g/3PV' _ g/3PV ^ 
or = 

Following [ij, we set A = 1 J- ?7 for infinitesimal p: 


Z[X^V,P,y]-Z[V,P,yl 

Z^[V,p,y]-Z[V,p,y]. 


( 7 ) 


ed^^DpPPV = [E, P, y] + dxZ^ [E, p, y] - Z[E, P, y] 


= dxZ^[V,P,y] 


A=1 


( 8 ) 


Z[V,P,y\(^J dr J d^x ^-0'1'VV - ^ J d^y pir.y) [{x - y) ■ X/sV{x - y)] p{t,x) 


where we’ve defined p{t,x) = 00t,£) 0(r,a;). Cancelling the partition functions on both sides, noting that thermal 
expectation values for the fields at the same r are independent of r so that the r integral pulls out a /3, and denoting 
the kinetic energy as KE: 


DPV = 2KE- J d°x J d°yp{T,y) [{x - y) ■ V sV{x - y)] p{t,x)^ 
which is the virial theorem in D dimensions (Eqs. (3.30) and (2.6) in Q and respectively). 


(9) 


^ In this paper we set h = m = 1. 
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III. N-BODY 

It is clear that this method can be generalized to the n-body case. Since by Eq. (2) the scaling transformation 
preserves f d^x (t, x)'ijj(T, x) (= f d^x p(t, x )), an n-body term transforms as 


1 

n! 


/ 



Xi p{t, Xi ) 





■ 


?;)• 


( 10 ) 


Setting V (xi,x^) = V{zcomj Z 2 , ■■■, Zn) where Zi = Xi — Xi and -?com is the center of mass of the x^’s gives 


ZQi O M ' zq o m ^ O M1 Z2 ■;••••) Zn ) 


DPV = 2KE- ^ y (^n ) 

For translationally-invariant systems, we can ignore the potential term in the 1st line. 


■ V?X(zCOM, 22, 

. 2=2 


( 11 ) 


IV. DIAGRAMMATIC PROOF OF P=P’ 

To prove diagramatically that the pressure P' corresponding to Z[X^V, /3, p\ is equal to the pressure P corresponding 
to Z[V,f3,p], it suffices to show that V,[X^V, f3, p] = X^fl[\^,p], where O is the grand potential. By the cluster 
expansion, O is given by the sum of connected vacuum graphs [ 3 . Using the Feynman rules, 0[1/, /3, p] oc {0)A4{f3, p), 
where <5^(0) expresses conservation of momentum of the vacuum and A4{l3,p) is the Feynman amplitude^ which is 
independent of U, since M. contains expressions like ^" 1 — f which in the continuum limit ( 2 -^)’^ f 

Taking (5^(0) oc U, it’s clear that 0[U,/3,^] oc VM{l3,p), so fl[X^V, P, p] = X^i}[V,P,p] in the continuum limit. 

Alternatively since Z[X^V, pi, p] = Z^[V, (3, p], another way to show P' = P is to show that the grand potential 
U''[U, /3, p] of Z^\y, j5, p] is larger by a factor of X^ than U[U, j5, p]. Then /3, p] = U[A^U, (3, p] = A^U[U, j3, p]. 

The grand potential is given by: 


= -l3-^\nZ\V,l3,p\. ( 12 ) 

By the cluster expansion, is given by the sum of connected vacuum graphs. Z^\V,f3,p] and Z\V,l3,p] have the 
same macroscopic parameters and only differ in that .Z^’s propagator is 




1 

"! 13 5 

XU)n 2A^ ^ 


(13) 


and that the potential is 


V^{x -y) = V (A {x -y)) 


(14) 


^ is the T-matrix, and 5^(0) = f ^-zO*x ^ y 

For finite volume, momenta are discrete and summed over: ki = Ani...An^ is a box of unit volume surrounding the discrete 

lattice point n*. In the limit of large L, / ^^ is assumed not to vary much, so any point within Ar7.i...An^ not on the lattice would 
still contribute the same value of / . Then = E ->■ / = / (f^/ (U)- 
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instead oiV(x—y). Fourier transforming Eq. (14) gives the relationship: 

( 15 ) 

The Feynman rules for the theory say that each vertex contributes its Fourier transform , where k is the 

momentum flowing through the vertex, and each propagator contributes Eq. (13). For vacuum graphs, all momenta 
k in the vertices and propagators are integrated over in loop momenta f Let us make the change of variables 

/ ( 27 r)^ “ / ( 271 -)^ relabel k as k. This will cause Ax{iuj,k) = A A{ioj,k) and 

I kj = in the loop integrals. 

Therefore, fl'*' is the same as fl, except for an overall scale factor of (A’^)^; where i/ is the number of vertices 

and L is the number of loops. Topologically, for connected vacuum graphs of the 2-body potential, L = ix + 1. So the 
overall scale factor becomes A^. Hence and therefore P' = P. 

This generalizes to translationally-invariant n-body potentials, and for spontaneous symmetry breaking. Suppose the 
interaction is of the form: 




(16) 


where m(i) is the number of fields in the interaction with spatial coordinate Xi, and M 
invariant potentials 


n 

^ m(z). For translationally- 
2=1 


X ^(!) 

A-DI"-!) ■ 


(17) 


So 




f I \ 

A^("-i) J 


(A^)^H. 


(18) 


Since L = — l) v + 1,® this again gives: 


= X^n. (19) 

For a diagram with a mixture of vertices of different types, T = ^ Vi + 1, where vi is the number of vertices 

i 

of type i, and Mi is the number of lines coming out of each vertex: 




n 

, i 

x^n. 


\Dr 


^£i£i \D{ni-l) 



Ai'i+i 

^ n 


( 20 ) 
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M lines come out of each vertex, and each line coming out is 1/2 of an internal line, so = I where I is the number of internal lines. 
The number of loops is the number of independent momenta, L = I — v 1. So L = — 1^ + 1. 
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V. SCALE EQUATION 

The virial equation, Eq. (9), can be recast into a different form that illustrates the effect of microscopic scales on the 
thermodynamics of a system. A simple way to see this is to write the potential as®: 


V{\x-y\) = 




( 21 ) 


\x-y\ 


/ is a dimensionless function whose arguments are the ratios of the couplings gi of V to their length dimension [gi] 
expressed in units oi\x — y \ (~ ig-g-p provides units of energy)^. Denoting r = \x — y\ 


r^ = -2U(r) + i- 

dr r dr 




= -2V{r) - -^[g^]g, 


(;&■) 


dgi 


= -2U(r)-[5,b|^. 


where the chain rule was used in line 2. Substituting this into Eq. (9) gives 


DPV = 2KE+ 2V - J d^x J d°yp{T,y) p{t,x) 

= p{t,x) 


( 22 ) 


(23) 


Rearranging: 


2E - DPV = 




(24) 


On the LHS of Eq. (24) are macroscopic thermodynamic variables. The RHS is a measure of the microscopic physics 
of the system. In particular, if the potential has no scales [gi] = 0 and no anomalies (i.e., T = IJ, you get 0 on the 
RHS, and Eq. (24) reduces to the equation of state for a non-relativistic scale-invariant system Q. 


VI. CONCLUSION AND COMMENTS 

The goal of this paper has been to highlight certain features in the derivation of the virial theorem for non-relativistic 
systems, which display a potentially important omission due to the presence of the Jacobian needed in the path- 
integral derivation developed here. Indeed, while we set J = 1 at the outset in order to make contact with the 
literature (specifically, Toyoda’s et al. work @-13)) Eq. (6) shows that the natural procedure would be to not assume 
this and keep the contribution of the Jacobian, regardless of whether or not there is a classical scaling symmetry. 
Obviously, in the latter case, one has to keep the Jacobian in order to incorporate the quantum anomaly as was shown 
in [J. The formal mathematical steps in the general case presented here are the same as in that paper, and Eq. (24) 
would become 


® We are now restricting ourselves to radial potentials. 

^ As an example, consider V{\x — y\) = — y\'^ + X\x — y\, where the coupling k has length dimension -4 and A has length dimension 

-3. Then / ^ ^|[fc] ; “ 1 |g—y ^ |g— couplings k and A provide the characteristic length scales. 
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2E - DPV = - 



where I 2 = 

matrix and functional trace). 



1 + X • V), and we have also used the 2x2 matrix notation of Q (Tr includes both a 


As with the work in [l| and Q, the key to assess the importance of the 
compute its contribution in detail, which implies a careful regularization 
tion. The actual details will depend of the type of potentials considered, 
generalization of these ideas. Work on this is currently in progress [^. 


Jacobian term rests upon one’s ability to 
procedure, and possibly also renormaliza- 
An interesting direction is the relativistic 
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